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Abstract 

The excitation spectrum of the one-dimensional spin-orbital model in a mag- 
netic field is studied, using a recently developed dynamical density matrix 
renormalization group technique. The method is employed on chains with up 
to 80 sites, and examined for test cases such as the spin-1/2 antiferromag- 
netic Heisenberg chain, where the excitation spectrum is known exactly from 
the Bethe Ansatz. In the spin-orbital chain, the characteristic dynamical re- 
sponse depends strongly on the model parameters and the applied magnetic 
field. The coupling between the spin and orbital degrees of freedom is found 
to influence the incommensuration at finite magnetizations. In the regions 
of the phase diagram with only massive spin and orbital excitations, a finite 
field is required to overcome the spin gap. An incommensurate orbital mode 
is found to become massless in this partially spin-polarized regime, indicating 
a strong coupling between the two degrees of freedom. In the critical region 
with three elementary gapless excitations, a prominent particle-hole excita- 
tion is observed at higher energies, promoted by the biquadratic term in the 
model Hamiltonian of the spin-orbital chain. 
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I. INTRODUCTION 



There is a wide variety of quasi-one-dimensional spin systems with short-ranged and 
quasi-long-ranged quantum liquid states. Examples for short-ranged compounds include 
the spin-Peierls chain, even-leg spin-1/2 Heisenberg ladders, and integer-spin Heisenberg 
chains. In contrast, half-odd-integer-spin chains and odd-leg spin-1/2 Heisenberg ladders are 
known to have a quantum critical ground state with quasi-long-range-ordered spin correla- 
tions. Many quasi-one-dimensional models, such as the antiferromagnetic spin-1/2 Heisen- 
berg chain with frustrating longer-range exchange interactions, have a phase diagram which 
contains both types of disordered ground state, depending on the parameters. Particular 
recent attention has focussed on spin-1/2 Heisenberg chains coupled to orbital degrees of 
freedom. |l]-|9| These systems are known to have a rich phase diagram, containing short- 
ranged and quasi- long-ranged regions, including an SU(4)-symmetric integrable point at a 
particular combination of the exchange parameters. 

In this work, we study the spin and orbital excitation spectra of spin-orbital chains in 
a magnetic field. By applying a newly developed dynamical density matrix renormalization 
group (DMRG) technique, [|K],[iT|] it is found that the admixture of the orbital degrees of 
freedom has a profound influence on both the field-dependent incommensuration, and the 
dispersion and spectral weight of the spin triplet modes, although only the spins couple 
directly to the applied magnetic field. Furthermore, this study is intended to provide useful 
information for inelastic neutron scattering studies on orbitally degenerate materials such as 
Na2Ti 2 Sb20 [TJ], where the d-electron of the Ti 3+ ion occupies an e g doublet, and NaX^Os, 
where orbital degeneracy appears to arise from the coupling of two types of VO5 chains. 
|l|JT3|-[l6| As long as experimentally available magnetic fields (~ 30 Tesla) are sufficiently 
large to overcome the respective spin gaps, our results for the dynamical properties of the 
incommensurate phases in the one- dimensional spin-orbital model may be tested on these 
or suitably related compounds. 

In this paper, the Hamiltonian of the one-dimensional spin-orbital model in a magnetic 
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field is studied. This can be written as 

N N N N 

H = Ji £ Si ■ S i+1 + J 2 ]T T, ■ T m + K£ (Si • S m ) (Ti ■ T i+1 ) - h £ S?, (1) 

j=l i=l i=l i=l 

where the spin-1/2 operators Si and the pseudo-spin-1/2 operators Tj describe the spins 
and orbits at a site i. |T7] The first two terms are the Hamiltonians for (pseudo-) spin- 1/2 



Heisenberg chains with exchange constants J\ for the spins and J 2 for the orbitals. The third 
term couples these two sets of degrees of freedom with a constant K. The fourth term is the 
Zeeman coupling of the spins to an external magnetic field aligned with the ^-direction. In 
the absence of a magnetic field, this Hamiltonian has an SU(4) symmetry for the particular 
parameter combination J\ = J 2 = K/A. This symmetry reduces to SU(2)xU(l) in an 
intermediate regime (0 < h < h c ), and to SU(2) at very high magnetic fields where the 
spins are completely aligned. If the parameter set ( J±, J 2 , K) is chosen away from the SU(4) 
point, the Hamiltonian remains invariant under SU(2)xSU(2). Furthermore, there is a Z 2 
symmetry about J\ = J 2 for vanishing magnetic fields. 

Recent studies have concentrated on the phase diagram and low-energy excitations 
Pjn| of H. At the integrable SU(4) point, the spectrum of excitations has been obtained 



from the Bethe Ansatz equations Effective low-energy theories were obtained from 



bosonization |||J, and approximate solutions from mean- field theory ||19|| . In addition, a 
variety of numerical techniques, including series expansions JE], numerical diagonalization 
||, Monte Carlo 0], and density matrix renormalization group 0, has been employed to 
analyze the ground state and low-energy properties of this Hamiltonian. 

In this work, a recently developed density matrix renormalization group technique is 
applied to study the full spin and orbital dynamics of the one- dimensional spin-orbital model 
in a magnetic field for arbitrary parameters (Ji, J 2 , K, h). The advantages of this method 
are that (i) relatively large Hilbert spaces can be studied (typical system size ~ 80 sites), (ii) 
the full excitation spectra can be obtained, and (iii) there are no constraints concerning low 
energies, small deviations from the SU(4) point, or the magnitude of the applied magnetic 
field. However, the spectra calculated with this method do exhibit some finite-size effects, as 



expected for low- dimensional quantum systems. Here we concentrate on the characteristic 
features of the dynamical spin and orbital structure factors for the various regimes of H. 
These quantities are defined as 

S(k,«) = -I toM» |^_J__ S .|» >, (2) 

r(*, u) = -iiimi m ^ |r ^ + . £ _ 1 g + go ^|^), (3) 

where j^o) is the ground state, E its energy, S% and are the Fourier transforms of the 
local spin and orbit operators S* and T*, and a small but finite broadening factor e = 0.05 Ji 
is chosen to represent the poles in S% and T%. 

In general, for the properties of the orbitally degenerate Hamiltonian H to be applicable 
to candidate inorganic materials such as Na2Ti2Sb20 and NaX^Os, the orbital energy lev- 
els must be nearly degenerate on the scale of the spin dynamics, which have characteristic 
energies in the range of J\ ~ 10-100 meV. While external magnetic fields couple directly 
to the spin degrees of freedom, applied external pressure could be a "dual" perturbation, 
primarily affecting the orbitals. It would therefore be of great interest to explore the dy- 
namical inelastic neutron scattering response in candidate spin-orbital compounds subject 
to a magnetic field as well as to external pressure. 

In the following section, we discuss briefly the numerical renormalization group method 
which is used to obtain the spin dynamics of the spin-orbital chain in a magnetic field. 
In section III, the field-dependent spectra are studied in the different parameter regimes 
(Ji, J2, K,h). A summary and conclusions are given in section IV. 

II. CALCULATION OF DYNAMICAL PROPERTIES WITH THE DENSITY 
MATRIX RENORMALIZATION GROUP TECHNIQUE 

The density matrix renormalization group (DMRG) method is known to be a precise 
numerical technique to calculate the ground state properties of quasi-one-dimensional sys- 
tems such as chains and ladders. pOf Recently, Hallberg ]10| and Kiihner and White JTT 



have proposed extensions of this method to extract dynamical correlation functions, based 
on continued-fraction expansions of the dynamical Green function. In this paper, we use a 
variant of this technique to calculate the momentum and frequency-dependent spectrum of 
the spin-orbital model in one spatial dimension. 

The zero-temperature, time-dependent correlation function of a momentum-dependent 
operator Ak is defined as 

A(k,t-t') = (y Q \AU(t)A k (t')\y ), (4) 

where |\I/o) is the ground state. In the following, A^ corresponds either to the spin operator 
Sf. or to the orbital operator T| . The Fourier transform of this correlation function can be 
expanded in the form of a continued fraction, 

AftH—ilm <*°l^l*°> , (5) 
71 uj + ie — a r- — 1 

where e — > is a small broadening factor. The coefficients a n and b n are then obtained from 
the recursion relation 

\f n+1 )=H\f n )-a n \f n )-b 2 n \f n - 1 ), (6) 

with |/o) = Afe|^o) an d b = 0. This procedure is very similar to numerical diagonalization 
calculations. However, there are specific technical problems intrinsic to the DMRG method 
which need to be addressed. pTi| . fTT| 

One complication arises when open boundary conditions (OBC) are used, which is usual 
in DMRG calculations to yield the highest possible numerical precision. In this case, the 
momentum is not a well-defined quantum number, and each momentum-dependent operator 
Ak generates wave packets of finite width, centered around k. To remedy this problem, Hall- 
berg has used periodic boundary conditions (PBC) to obtain a basis with a good momentum 
quantum number, whereas Kiihner and White have used OBCs and the Parzen filter to elim- 
inate the artificial surface excitations on the surface boundaries. In this paper, we adapt 
OBCs combined with the Parzen filter to obtain the full dynamical spectra associated with 
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S% and T%. OBCs are known to give more precise ground states than PBCs, and therefore 
larger clusters can be studied, while the width of the corresponding wave packets is greatly 
reduced. 

Further improvement of the dynamical DMRG method may be obtained by using more 
target states besides the ground state, implemented with "correction vectors". [TT|J2T[ We 
find that a rather large cut-off dimension of the density matrix is needed to reach convergence 
in the spin-orbital chain for the target states at higher energy. These convergence problems 
tend to be exacerbated for excitations away from the dominant wave vectors. However, 
with the algorithm used in this work, the full structure of the excitation spectrum can be 
obtained with rather high accuracy, even if only the ground state is kept as a target state, 
although the relative overlap (matrix elements) with higher-energy states is reduced. 

We have tested the algorithm on the spin-1/2 antiferromagnetic Heisenberg chain in 
a magnetic field. The Hamiltonian for this system corresponds to the parameter choice 
( Ji = J, J 2 = K = 0) for the spin-orbital model. One observes that the spin-triplet excitation 
spectrum of this simple model matches very well the exact results from the Bethe Ansatz. 
P^] A more detailed discussion of this spectrum will be presented in the next section. Chains 
with N = 80 spins were studied. The cut-off factor, i. e. the number of states kept per block, 
was taken to be m = 160. We have also used the Lanczos method on lattices of up to 16 
sites to check the accuracy of the DMRG results for other non-trivial parameter sets of the 
spin-orbital chain. 



III. DYNAMICAL SPIN STRUCTURE FACTOR OF THE SPIN-ORBITAL 

CHAIN 

The rich phase diagram of the one-dimensional spin-orbital model, represented by the 
Hamiltonian H in Eq. 1, has been the subject of a number of recent studies. |1|-|3| The 
main result of these works is summarized in Fig. 1, which we have adapted from Ref. 
||. The various phases can be characterized by spin and orbital quantum numbers (S, S z ) 
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and (T, T z ), and by effective Luttinger parameters in each channel, such as gaps (As, Ao), 
soliton velocities (v$, vo), and radii of compactification. For finite values of the spin-orbital 
interaction K , the spin and orbital degrees of freedom are coupled, and the excitation spectra 
generally contain strongly admixed spin and orbital modes. 

Let us first briefly describe the phase diagram of H . Phase I is "ferromagnetic" in 
the spin and orbital channels, i.e. both degrees of freedom are fully polarized, and there 
is an infinite degeneracy in the ground state. In phase II, the spin degrees of freedom 
remain polarized, whereas the orbitals are in a quasi-long-range-ordered antiferromagnetic 
state with a characteristic orbiton velocity Vq- In phase III, the roles of spins and orbitals 
are interchanged. Phase IV is fully gapped in all channels (As 7^ 0, A G 7^ 0) because of 
spontaneous dimerization. In particular, this phase contains a Majumdar- Ghosh point at 
J\ — J2 — 3K/ 4 which has a matrix-product ground state. [ 23 124] It has recently been shown 
that phase V contains 3 gapless excitations with non-universal, parameter-dependent 
soliton velocities. At the SU(4)-symmetric point (filled circle in Fig. 1), these three modes 
have one common velocity, v = nJi/2 = 7rJ 2 /2 = ttK/8. |2B[] 

Let us now turn to the excitation spectra of H in these regimes, starting with phase III 
where the orbital degrees of freedom are fully polarized. In Fig. 2, the spectrum traced by 
the dynamical spin structure factor is shown on a lattice of 80 sites, using the parameter 
set (Ji = — J 2 = K/2). The poles have been given a width of e = 0.05Ji, and the color 
scheme highlights the areas of high spectral intensity with bright colors. Since in this phase, 
the orbital degrees of freedom couple only trivially with the spins, the analysis turns out to 
be particularly simple. The spin excitation spectrum of this system corresponds to a spin- 
1/2, antiferromagnetic Heisenberg chain with exchange coupling Ji, which is dominated by 
a low-energy 2-spinon continuum (see Fig. 2(a)). |27| The elementary excitations of this 
system are spin-1/2 solitons with dispersion u a (q) = 3nJi sin (g)/4 for q G [0, 7r]. Thus the 
2-spinon continuum has a lower cut-off at u>i(k) = 37rJi| sin (A;)|/4 and an upper cut-off at 
w u {k) — 37rJi| sin (k/2)\/2. P2 , |28" | The dominant singularity of the dynamical spin structure 
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factor at wave vector ir diverges according to S(tt,u) oc Jin {oS)/u as lu — > 0. 

In a finite magnetic field, this mode becomes incommensurate with nesting vectors 2ftF = 
7r(l ± 2m), where m G [0,1/2] is the magnetization. This behavior is confirmed by our 
dynamical DMRG calculation, as shown in Fig. 2(b). Here, a finite magnetization, m = 
4/40, leads to an incommensuration with 2kp = 127r/10 and 87r/10. Because of the relatively 
large system size (N = 80), the spectra appear to be almost continuous, as expected in the 
thermodynamic limit. However, finite-size effects are still visible, causing gaps and uneven 
spectral weight distributions at higher energies. Furthermore, the momentum is quantized 
in 80 slices of size 7r/40. In spite of these limitations, the agreement of the numerical results 
with the exactly known spectra is encouraging, and gives confidence to proceed to parameter 
regimes with more complex excitation spectra. 

At the SU(4) point (J± = J2 = K/4), the spectrum is also known exactly from the 
Bethe Ansatz. A comparison of the spectra obtained using the DMRG method with the 
exact results can thus help to further evaluate the accuracy of the numerical approach. In 
Fig. 3, the spectrum traced by the dynamical spin structure factor is shown on a lattice of 
80 sites. From the solution of the Bethe Ansatz equations, it is known that the elementary 
excitations of the SU(4) chain are 4-fold degenerate spin- 1/2 spinons (a) and pseudo-spin- 1/2 
orbitons (r), and 6-fold degenerate spin-1 or orbital-1 solitons (u). [18|]. The corresponding 
dispersions of these three elementary excitations are 

Uff (q) = [V2 cos {-q + 3tt/4) + l] , q 6 [0, 3vr/2] 
Jn 

w T (?) = -y [V2 cos (-q + tt/4) - l] , q e [0, tt/2] (7) 
Jtt 

u v(q) = -y cos (-q + it/2), q e [0, tt] 

where J = J\ = J2 = K/4. The spin-triplet excitation spectrum of an SU(4) chain with 
N = 4n sites, shown in Fig. 3(a), contains convolutions of elementary a and r excitations, 
as well as v excitations. The z-component of the dynamical spin structure factor S z (k,u) 
couples directly to the 15-fold multiplet of cr — r pairs, the 20- fold multiplet of 2v excitations, 
and with much smaller matrix elements to the 45-fold 2a — v and the 35-fold 4a multiplets. 
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Fig. 3(b) illustrates that the spectrum obtained with the dynamical DMRG traces 
exactly the two-soliton continua, spanned by uj a (q) + oj T (q') and by cj u (q) + LU v (q'). These 
continua have characteristic low-frequency power-law singularities (S z (k, u) oc uj~ a ) at the 
Fermi vectors tt/2 and 37r/2, and a weaker divergence at tt. [21]] The relatively high spectral 



intensity seen at larger frequencies and centered around the wave vector tt, can be attributed 
to a particle-hole excitation in the z/-channel which will be discussed later. We note that the 
accuracy of the dynamical DMRG method is known to deteriorate at larger energies, and 
hence the spectral intensities are less reliable in this regime. 

When an external magnetic field is applied along the ^-direction, the spectrum of the 
SU(4) chain also becomes incommensurate. |J This effect can be understood quantitatively 
by considering the splitting of spin-up and spin-down bands due to the applied field. The 
corresponding nesting vectors depend on the magnetization m e [0, 1/2] as 2kp^ = 7i(l + 
2m)/2, 2fcp| = 7r(l — 2m)/2, and 4/cp-f = 4fcp^ = 7r(l — 2m). This low-frequency behavior is 
reflected in the finite-size DMRG spectra, shown in Fig. 4, with magnetizations m = 4/40 
and m = 12/40. For example, with m = 4/40, the nesting vectors are 67r/10, 47r/10, and 
87r/10, as shown in Figs. 4(a) and (c), whereas for m = 12/40 they are 87r/10, 27r/10, and 
47r/10 respectively (Figs. 4(b) and (d)). 

Furthermore, the overall spectral weight is reduced as the phase space of spin-triplet 
excitations shrinks with increasing field. At finite magnetic fields, there remains a commen- 
surate soft mode at wave vector it in the orbital z/-channel which does not couple to S(k, uo). 
It is therefore interesting to examine the orbital excitation spectrum traced out by T(k,u). 
The orbital dynamical structure factor of the SU(4) chain in a magnetic field is shown in 
Figs. 4 (c) and (d). In the absence of a magnetic field, the orbital spectrum is identical to 
the spin spectrum (Fig. 3). In a finite magnetic field there is one orbital mode, the spin-1 
z/-channel, which does not become incommensurate, but stays soft at wave vector ix for all 
magnetic fields. On the other hand, the strongly mixed a — r pairs show incommensuration 
with the same nesting vectors as for the spin channel. The dominant a — t singularities at 
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low frequencies occur at 2fcp-|- for the orbitons and at 2kp[ for the spinons. 

The excitation spectra at the Majumdar- Ghosh point (Ji = J 2 = 3i^/4) are shown in 
Fig. 5. The zero-field ground state at this point in the dimerized phase IV is a doubly 
degenerate product of spin and orbital singlets with an energy of — 3Ji/4 per site and a 
gap of approximately 3K/8 to the lowest excitations. []I|, [2"3T , p4| The spectrum contains a 
2-soliton continuum of states with minima at wave vectors k = and tt, and maximum 
spectral weight at wave vector tt. |24[] In addition, there are massive magnon excitations 



corresponding to Haldane triplet bound states. [2£| These are centered at k = tt/2 and 



37r/2, indicating a doubling of the real-space unit cell. In Fig. 5(a), we have eliminated the 
artificial low-energy states which arise from surface excitations at the open boundaries of 
the 80-site cluster. These states have vanishing spectral weight in the thermodynamic limit. 

A finite magnetic field leads to the deconfinement of bound spinons. |3(J This is man- 
ifested in a lowering and eventual vanishing of the corresponding spin gap when the field 
exceeds a critical value, which for the Majumdar- Ghosh point is at /i cl pa 3K/8. In Fig. 
5(b), the spin excitation spectrum is shown in this incommensurate regime (h c i < h < /i c2 ) 
at a magnetization m = 4/40. The dominant low-energy singularities of the spin-1 magnons 
occur at wave vectors 4kp = 7r(l ±2m). The 2-spinon continuum at higher energies appears 
to be almost unaffected by the finite magnetic field. Interestingly, a strong low-energy singu- 
larity at k = tt in the orbital channel is induced by the applied field, as observed in Fig. 5(c), 
which is indicative of an effective attractive potential between the incommensurate spin-1 
and the commensurate orbital- 1 z/-type solitons in this case. This magnetic-field- induced 
two-component Luttinger Liquid behavior, with the orbital channel remaining commen- 
surate and the spin channel becoming incommensurate, has also been observed in recent 
bosonization studies of the spin-orbital chain. 

Finally, let us turn to the gapless phase V. In Figs. 6 and 7, the spin and orbital 
excitation spectra are shown at J\ = J2 = K/8 and at J\ = J2 = 0. As K — > 0, a prominent 
high-energy excitation emerges, centered at wave vector tt and frequency uj pa 0.8K. This 
excitation carries the quantum numbers (S,T) = (0,1) and (1,0) which indicates a particle- 

10 



hole-type bound state, involving a spin singlet and an orbital triplet or vice versa. From 
finite-size scaling, we observe that the energy gap between the ground state and this level 
increases with system size. If J\ = J 2 = 0, the eigenstates of the K term can be found 
exactly on a 4-site cluster. The ground state with energy — 3K/A and the excitons with 
energy are in this case 
4 

\ GS ) = ^(-i)* I (^0^+2,0X^,^+2,5) - (Si,OS i +2,o)(tt,S t T+2,s) 
1=1 

+ (^+l,0^i+3,o) (^5^+2,5) ~~ ( t to t i+2,o)( t i+l,S t t+3,s)) 
4 4 

k S ) = £ l(^oV+2,o)(^5 s i+2,5)) and |tt, 0) = J2 I(^sV+2,s)(^o s *+2,o)> 
1=1 1=1 

where, for example, t/" 5 denotes an up-up triplet state in the spin channel, involving spins 

at sites i and i + and s i+2 ,o denotes a singlet in the orbital channel at sites % + 2 and i + 3. 

This high-energy bound state is observed in the entire phase V, but is most prominent in 

the strong- K limit. 

The magnetic-field dependence of the spectra in phase V can be analyzed analogously to 
the SU(4) limit. In particular, the zero-energy modes can be understood within the band 
picture discussed above, with nesting instabilities at 2k.pi, an< ^ 4/cpf=4A;^. The finite- 
size gaps in these spectra appear to be most obvious at the critical point, and become less 
severe as Ji, J 2 — >• 0. 

IV. SUMMARY 

In summary, the spin and orbital excitation spectra of the spin-orbital chain were studied 
in various regimes of the model Hamiltonian H, given in Eq. 1. The characteristic dynamical 
response depends strongly on the parameter set (Ji, J 2 ,K) and on the applied field h. In 
phases II and III with either fully polarized spin or orbital degrees of freedom, the remaining 
gapless mode has the characteristic spectral response of a (pseudo)-spin-l/2 Heisenberg 
chain. On the other hand, in the dimerized regime (phase IV), a finite magnetic field is 
needed to close the spin gap. Interestingly, an additional orbital mode is found to become 
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gapless at the commensurate wave vector 7r in the partially spin-polarized regime (h c i < 
h < h c2 ), whereas the Fermi momenta of the incommensurate gapless spin modes depend 
on the magnetization via 2k f = 7r(l — 2m). In the gapless regime (phase V), there are three 
elementary gapless excitations. The dynamical response consists of various overlapping two- 
particle continua, arising from pairs of these elementary excitations. At higher energies, a 
particle-hole bound state whose origin lies in the biquadratic i^-term, is found to dominate 
the spectrum. 

We are greatful to Andreas Honecker, Till Kiihner, Bruce Normand, and Edmond Orignac 
for useful discussions, and we thank Till Kiihner for valuable help in developing the dynam- 
ical DMRG algorithm. S. H. acknowledges the Zumberge foundation for financial support. 
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FIG. 1. Phase diagram of the one-dimensional spin-orbital model, adapted from Ref. 3. Phase 
I: spins and orbitals are fully polarized. Phases II and III: one of the channels is fully polarized, 
whereas the other is a quasi-long-range-ordered antiferromagnet. Phase IV: dimerized regime. 
Phase V: gapless regime. Circle: SU(4)-symmetric point. 
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FIG. 2. Spin excitation spectrum of the spin-orbital chain at (Ji = — J2 = K/2), obtained 
with dynamical DMRG on a 80-site cluster, (a) zero magnetic field m = 0, and (b) m = 4/40. 
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FIG. 3. (a) Dynamical spin structure factor of the SU(4)-symmetric spin-orbital chain. 
This spectrum was calculated using the dynamical DMRG technique on an 80-site cluster, (b) 
Spin-triplet spectra composed of a — r pairs (dots) and 2v pairs (pluses), obtained from convolu- 
tions of the single-soliton dispersions, which are known exactly from the Bethe Ansatz solution. 
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FIG. 4. Spin and orbital excitation spectra of the SU(4) chain in a magnetic field, (a) (upper 
left) spin dynamics at m = 4/40, (b) (upper right) spin dynamics at m = 12/40, (c) (lower left) 
orbital dynamics at m = 4/40, and (d) (lower right) orbital dynamics at m = 12/40. 
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FIG. 5. Spin and orbital excitation spectra of the spin-orbital chain at the Majumdar- Ghosh 
point (Ji = J2 = 3if/4). (a) (upper) spin and orbital dynamics at m = 0, (b) (lower left) spin 
dynamics at m = 4/40, and (c) (lower right) orbital dynamics at m = 4/40. 
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FIG. 6. Spin and orbital excitation spectra of the spin-orbital chain at {J\ = J2 = K/8). (a) 
(upper) spin and orbital dynamics at m = 0, (b) (lower left) spin dynamics at m = 4/40, and (c) 
(lower right) orbital dynamics at m = 4/40. 
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FIG. 7. Spin and orbital excitation spectra of the spin-orbital chain at {J\ = J2 = 0). (a) 
(upper) spin and orbital dynamics at m = 0, (b) (lower left) spin dynamics at m = 4/40, and (c) 
(lower right) orbital dynamics at m = 4/40. 
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